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Abstract 

The strong similarity between the magnetohydrodynamic (MHD) turbulence and initially 
isotropic turbulence subject to rotation is noted. We then apply the MHD phenomenologies 
of Kraichnan and Matthaeus k Zhou to rotating turbulence. When the turbulence is subject 
to a strong rotation, the energy spectrum is found to scale as E(k) = Cu(9.eY /2 k ~ 2 , where 
n is the rotation rate, k is the wavenumber, and e is the dissipation rate. This spectral 
form is consistent with a recent letter by Zeman. However, here the constant Cq is found 
to be related to the Kolmogorov constant and is estimated in the range 1.22 - 1.87 for 
the typical values of the latter constant. A ‘rule’ that relates spectral transfer times to 
the eddy turnover time and the time scale for decay of the triple correlations is deduced. A 
hypothesis for the triple correlation decay rate leads to the spectral law which varies between 
the ‘-5/3’ (without rotation) and ‘-2’ laws (with strong rotation). For intermediate rotation 
rates, the spectrum varies according to the value of a dimensionless parameter that measures 
the strength of the rotation wavenumber k Q = (f^/c) 1 / 2 relative to the wavenumber k. An 
eddy viscosity is derived with an explicit dependence on the rotation rate. 
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1 Introduction 


For homogeneous, isotropic, statistically steady turbulence, the classical Kolmogorov phe- 
nomenology leads to the well known energy spectrum 


m = cxtWk-w. (1) 

Here e is the energy dissipation rate, k is the wavenumber, and C K is the Kolmogorov 
constant. The Kolmogorov hypotheses implies that the energy-containing range excitation 
does not affect energy transfer within the inertial range. Therefore, the average rate of energy 
dissipation is identified with the rate of spectral energy transfer and the rate of energy input. 
The nonlinear (or eddy turnover) time-scale, T n ,(k) = \k-FAk)}-\ is then equivalent to the 
spectral transfer time, t s (Kolmogorov 1 ; Batchelor 2 ; Monin and Yaglom 3 ). 

In this report, attention is called to the strong similarity between the magnetohydrody- 
namic (MHD) turbulence and isotropic turbulence subject to solid body rotation. For MHD 
turbulence, Kraichnan 4 argued that the propagation of the Alfvenic fluctuations disrupts 
phase relation and thereby may be expected, on the average, to decrease energy transfer. 
Experimental studies 5 " 8 , direct numerical simulations (DNS) 9 " 12 , large-eddy simulations 
(LES) 9,13 , and closure approximations 14 have established that solid body rotation suppresses 
the nonlinear energy cascade from large scale to small scale. Specifically, uniform rotation 
causes plane waves to propagate with phase speed 2 Mc z /k. Here the rotation vector is consid- 
ered to act along the vertical (ft). The effect of rotation on the spectral transfer is through 
phase scrambling. Based on this evidence, the phenomenology 4 . 15 originally developed for 
MHD turbulence is applied to rotating turbulence. 

2 Energy Spectra and Transfer 

The Kolmogorov’s phenomenology cannot be simply extended to more complex situations 
which may be governed by other time-scales. Special generalizations of the theory are nec- 
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essary. Kraichnan 4 generalized the Kolmogorov theory for MHD turbulence so that the 
influence of the Alfven wave on the MHD spectrum can be included. In addition to the 
nonlinear time scale, there is a much shorter Alfven time scale = [kB 0 }~\ where B 0 is 
the characteristic magnetic field strength in Alfven speed unit. In order to infer the form of 
the inertial-range spectrum, it is necessary to estimate the magnitude for the triple correla- 
tions. If the energy-containing range excitations were absent, it would be expected that the 
nonlinear interaction would build up substantial triple correlations in the local dynamical 
time r n i(k). With the energy-containing excitation present, however, the Alfven time scale 
Tji is the effective time for relaxation of the locally built-up phase correlations through prop- 
agation. In general, r 3 , the time scale for decay of triple correlations which is responsible 
for inducing turbulent spectral transfer, may depend on any relevant turbulence parameters. 
Because energy is conserved by the nonlinear interaction and a local cascade has been as- 
sumed, e is independent of k. Local cascade also implies that e is explicitly proportional to 
r 3 and depends on the wavenumber and on power of the omni-directional energy spectrum 
A simple dimensional analysis leads to 

e = A 2 T 3 {k)k 4 E 2 (k) (2) 

where A is a constant. In steady state isotropic turbulence, the only time scale available for 
r 3 (k) is the nonlinear time scale, r n/ ( *); therefore the Kolmogorov 1 spectrum is recovered. 
Using r 3 (fc) = r^(fc), Kraichnan 4 obtained the MHD energy spectrum 

E(k) oc ( B 0 e) 1,2 k- 3/ 2 . ( 3 ) 

Recently, Squires et al. 13 extended one point analysis of the long-time evolution of kinetic 
energy and length scales to rotating homogeneous turbulence. In regime of high Reynolds 
numbers and low Rossby numbers, rotating turbulence is characterized by two disparate time- 
scales: a short timescales associated the rotating frequency r fi = 1/D and a nonlinear time 
scale. Squires et al. 13 assumed that the correlation time of the nonlinear triadic interactions 
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is directly proportional to the short time scale r„. The predictions of the asymptotic decay 
for' the kinetic energy in rotating turbulence, based on this crucial assumption on the time 
scale, have been confirmed within a few percent by their LES runs for times of order O(10 3 ) 
of initial turbulence. A direct application of r 3 = r„ produces the following energy spectrum 
for turbulence subject to strong rotation: 


E(k) = Cn(ne)'' 2 k- 2 , (4) 

where C u is a constant. The spectral form in Eq. (4) is consistent with a recent letter of 
Zeman 16 which motived present study. However, here we will determine the constant C n . 

For steady isotropic turbulence, the wavenumber dependent spectral transfer time r s (ifc) 
can be defined through e = U 2 k /r s (k), where U k = [*£(*)JV* is the characteristic velocity of 
eddies at wavenumber k and e is the steady, wavenumber-independent energy transfer rate. 
Using (2) and the definition of r nl and U k , it was found that 


Ts (k) = IMS 

1 ’ A 2 Tn(k) ■ 


( 5 ) 


A corollary to this development is the time for the spectral energy transfer time is increased 
to a value greater than r„,; thus the nonlinear energy transfer is suppressed by the rotation. 
Similar suggestion has been made for MHD turbulence 17 " 18 . 

For sufficiently large Reynolds numbers for which a universal equilibrium exists, the 
expression for the skewness 


g. _ 3(30) 1 / 2 / 0 °° k 4 E(k) dk 

7 U [f™k*E{k)dk}W ( 6 ) 

is exact (Batchelor 2 , Monin and Yaglom 3 ). Eqs. (4) and (6) implies that the skewness 

factors of the spatial derivatives of the fields approach to zero as the rotation rate go to 

infinity (S 3 a 1/fi / ). From the single point modeling perspective 19 , the reduction of the 

skewness leads to the decreasing of the production term in the dissipation equation. 

Following the MHD development, we make use of the viewpoint that the lifetime of 

triple correlations in rotating turbulence might be more accurately treated by taking into 
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account the possibility that these correlations decay because of the influence of both wave 
propagation and nonlinear triadic interactions (Matthaeus and Zhou 15 ). The simple choice 


r 3 {k) T n i(k) + T Q (k) 

satisfies the appropriate limiting cases: r 3 (k) -*■ T ni without rotation and r 3 {k) -»• t q with 
strong rotation. The most general rule is 


r s {k) = 


1 

A 2 r 3 (k ) 


(8) 


Substituting (7) into (2), the generalized inertial range energy spectrum is 


E{k) = Z 2 A~ 4/3 e 2/3 k 5/3 , 


( 9 ) 


where Z is the solution of the dimensionless equation 


Z 4 = Z 0 + Z 


(10) 


where 

7 _ 42/3 g = r^ll2/3 (11) 

Z ° _ A (ek 2 ) 1 / 3 k J 

The parameter k u = (0 3 /e) 1/2 has appeared in Zeman 16 . Eqs. (9)-(ll) reduce to the non- 
rotating Kolmogorov “-5/3” spectrum when Z 0 -» 0 (so that Z -» 1), and to our rotating 
modified “-2” spectrum when Z 0 -► oo (so that Z -» Zq / 4 ). As shown in Matthaeus and 
Zhou 15 , Eqs. (10)-(11) have only one physically acceptable solution to this quartic equation: 

z = ^Vy + y-^ + 2 \/^ 2 + 4 ^ ) (12) 


where 


Y = 


N 2 


1 / 1 A%o 3 

2 + V4 + ( T-» + \ 


1 /l /4Zo .3 

i-Vi + 'T 1 ' 


(13) 


From Eqs. (9)-(ll), the constant of the rotating spectrum can be estimated. Taking the 
Kolmogorov limit, we find that A = C^ /A (Matthaeus and Zhou 15 ). The strong rotation 
limit then leads to C Q = 1/A = 1.22-1.87 for the typical range of Kolmogorov constant. Fig. 
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1 illustrates spectra obtained this way with e = 1 and A = 0.7 (we have chosen C K = 1.5), 
varying the reference rotation rate over values ft = 0, 10, and 100. The zero rotation case is 
a pure Kolmogorov spectrum while the ft = 100 case is very nearly k~ 2 spectrum (Eq.4). For 
intermediate rotation rates the spectrum varies smoothly between these two limiting forms, 
according to the increase of the controlling parameter Z 0 with increasing ratio k Q /k. 

3 Eddy Viscosity 

The rotation dependence eddy viscosity can be estimated using Zeman’s procedure 16 . Fol- 
lowing Tennekes and Lumley 20 , it is reasonable to assume that in the inertial range 

T(k ) = e oc v T (k")S 2 (k') ( 14 ) 

where, approximately, k / 3 < V < k and k" = 3 and S(k) = 1 /r nl = [&E{k)]W (Ref. 21 

verified the predictions of Tennekes and Lumley 26 on the scaling laws of the energy transfer) 
Using (2), the spectral eddy viscosity is 

= (15o) 

where is a constant and y = PE(k)/(l\ following Zeman« Equation (15a) is expected 
to find use in the spectral LES of rotating turbulence. 

The inertial range wavenumber k can be related* 2 -*" to the turbulent kinetic energy Ii 
and dissipation rate e 

k = f^( 3 ^/2) 3/2 - 

For the inertial range spectrum, Eq. (15a) can be rewritten in physical space: 


vt{x) v 0 ^ 1 + 0 36/a) / e * (156) 

The eddy viscosity form (15b) is only appropriate for homogeneous turbulence without mean 
velocity gradients. For turbulent flows with mean velocity gradients in a rotating frame, 
the effect of rotations appears in conjunction with the mean vorticity (Speziale 25 ) and Eq. 
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(15b) can be extended in a straightforward way by replacing ft by {}V l 3 W l 3 f 12 (Gatski and 
Speziale 26 ). Here W {j = U l3 + ae mjt ft m is the extended absolute vorticity tensor defined m 
Gatski and Speziale 26 (o is a constant). The constant v ' 0 can then be evaluated in homoge- 
neous shear flow where Wij reduces to the uniform shear rate. Such flows are deferred to a 

future study. 


4 Discussion 

These results are now compared to Zeman 16 . He started with a simplified Lagrangian de- 
scription of the relation between stress and mean strain in rotating turbulence and proposed 
the modified the inertial range spectrum when t u > r n( 

E(k) a e 2 ' 5 ft 4/5 r n / 5 . (16) 


Zeman’s generalized energy spectrum is obtained by solving 


1+2 y + y 2 


(17) 


where x = k/k Q . The rate of energy transfer is given by 

_ 3/2 kUf 

fc 1 + Cs/y- 


(18) 


Comparing (18) with (2), we found that Zeman’s corresponding expression of the time scale 
of the triple velocity correlation is 


?3 = C* 3 [ 1 + y] 1 7nl i 


(19) 


where C 3 is a constant introduced by Zeman. 

Under rapid rotation, Eq. (19) reduces to 

m 1 t A 

T3W C 3 C^r n/ (fc)' 


( 20 ) 


Eq. (20) indicates that Zeman’s time scale for the triad nonlinear interactions differs from the 
key assumption of Squire et al. that the correlation time of the nonlinear triadic interactions 
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is directly proportional to the short time scale r n = 1/(2. Therefore, it appears that Zeman ’s 
basic assumption that leads to his energy spectrum is not supported by LES computation. 
Indeed, Zeman 16 noted that a further assumption 

1 +C 4 y~ 1 =y~ 1/2 (21) 

is needed to obtain the “-2” spectrum (Eq. 4) from Eqs. (17)-(19) (C 4 is another constant 
introduced by Zeman). 

Bershadskii et al. 2 ‘ have proposed a E(k) oc £ 2/3 |fc| - '/ 3 spectral form for rotating turbu- 
lence with some observational supports from stratospheric turbulence. Here ( is the rate of 
spontaneous helicity generation. This spectral form was originally obtained by Brissand et 
a/. 28 for the case of three dimensional turbulence, Bershadskii et al. 27 argued that the quasi- 
two-dimensional turbulence regime arises as a result of a spontaneous breaking of reflexional 
symmetry, which in turn is a consequence of the instability of two dimensional turbulence 
to three dimensional helical traveling waves and solitons. 

5 Conclusions 

The strong similarity between the MHD turbulence and isotropic turbulence subject to rota- 
tion has been noted. The MHD phenomenologies of Kraichnan 4 and Matthaeus and Zhou 15 
were applied to rotating turbulence. When the turbulence is subject to strong rotation, an 
energy spectrum, E{k) = C Q (neY' 2 k~ 2 , where the constant C Q = 1.22-1.87 is found. A hy- 
pothesis on the triple correlation decay rate, Eq. (7), leads to the spectral law given in Eqs. 
(9)-(ll), which varies between the ‘-5/3’ (without rotation) and ‘-2’ laws (with strong ro- 
tation). We deduced a ‘rule’ (Eq. 8) that relates spectral transfer times to the eddy turnover 
time and the time scale for decay of the triple correlations. For intermediate rotation rate, 
the spectrum varies according to the value of a dimensionless parameter that measures the 
strength of the rotation wavenumber = (fP/e) 1 / 2 relative to the wavenumber k. An eddy 
viscosity is derived which explicitly depends on the rotation rate. 
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